Entanglement Mean Field Theory and the Curie- Weiss Law 
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The mean field theory, in its different hues, form one of the most useful tools for calculating 
the single-body physical properties of a many-body system. It provides important information, 
like critical exponents, of the systems that do not yield to an exact analytical treatment. Here we 
propose an entanglement mean field theory (EMFT) to obtain the behavior of the two-body physical 
properties of such systems. We apply this theory to predict the phases in paradigmatic strongly 
correlated systems, viz. the transverse anisotropic XY, the transverse XX, and the Heisenberg 
models. We find the critical exponents of different physical quantities in the EMFT limit, and in the 
case of the Heisenberg model, we obtain the Curie- Weiss law for correlations. While the exemplary 
models have all been chosen to be quantum ones, classical many-body models also render themselves 
to such a treatment, at the level of correlations. 



I. INTRODUCTION 

Characterizing many-body systems by understanding 
their different phases is a key issue in physics. However, 
it is only in a few cases that exact analytical techniques 
can be applied pL] . It is therefore crucial to have approx- 
imate methods to deal with such systems to predict their 
different physical properties [2|. A very useful method 
is to use mean field theory (MFT) which ren- 
ders the many-body physical system into one with a sin- 
gle particle. MFT allows one to predict the single-body 
physical properties, like magnetization, susceptibility, of 
the system. Its importance lies in the fact that these 
MFT-reduced one-body properties can correctly predict 
the thermal fluctuation- and quantum fluctuation- driven 
phase transitions of the system, as well as the critical ex- 
ponents of the such single-body physical quantities [1, Q . 

Useful as it is, there are important limitations of a 
mean field theory a first being that it is not possi- 
ble to predict the multi-party physical properties, like 
entanglement, of the system by using this theory. 

In this paper, we propose an "entanglement mean field 
theory" (EMFT), which transforms an interacting many- 
body physical system into a two-body one, while still 
retaining certain footprints of the interactions in the 
many-body parent, using which it is possible to calcu- 
late the two-body properties, like two-point entanglement 
and two-point correlations, of the system, and predict 
the critical phenomena in it. Moreover, it is possible to 
calculate the critical exponents of two-body properties. 
The theory can be applied for detecting phase transi- 
tions driven by thermal fluctuations as well as quantum 
fluctuations. At the same time, both quantum as well 
as classical models can be treated by EMFT, where in 
the classical case, this will be only until the level of cor- 
relations. In this paper, we will only consider the ap- 
plications of EMFT to quantum systems. We will con- 
sider three paradigmatic classes of interacting spin mod- 
els, viz. the transverse anisotropic quantum XY (which 
includes the transverse quantum Ising), the transverse 
quantum XX, and the quantum Heisenberg models. The 
phase transitions of these models, both temperature- 



induced phase transitions as well as zero temperature 
quantum phase transitions, are faithfully signaled by the 
corresponding EMFT entanglements. The Curie- Weiss 
law of susceptibility is an important prediction of the 
mean field-reduced Heisenberg model in the paramag- 
netic regime Q. We indicate the corresponding law for 
EMFT-reduced correlations by using the EMFT-reduced 
Heisenberg model. 

The entanglement mean field theory opens up the pos- 
sibility of investigating the behavior of entanglement and 
other two-body properties of many-body systems, par- 
ticularly for the ones which does not lend themselves to 
an analytical treatment. At the fundamental level, this 
forms, potentially, an important link between many-body 
physics and quantum information science Q. 



II. THE ENTANGLEMENT MEAN FIELD 
THEORY: XY MODEL 

Before presenting the entanglement mean field the- 
ory, let us briefly remind ourselves the mean field the- 
ory. Consider the transverse quantum anisotropic XY 
model with nearest-neighbor interactions, described by 
the Hamiltonian 



(1) 



which represents a system of interacting spin-1/2 par- 
ticles on a d-dimensional cubic lattice. The coupling 
strength J/2 is positive, the anisotropy 7 S (0,1], and 
the transverse field strength h is also positive, ax, cry, 
and CTz are the Pauli matrices for the spin degree of free- 
dom of a spin-1/2 particle, (ij) indicates that the cor- 
responding sum runs over nearest neighbor lattice sites 
only. The mean field theory consists in assuming that a 
particular spin, say iQ, is special, and replacing all other 
spin operators by their mean values. Denoting the mean 
values of the spin operators ax, ay, a^ as nix, 'my, Wz, 
respectively, this leads to an MFT Hamiltonian H^j^j, 
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0,3 ■ We then solve the self-consistency equations (mean 
field equations) 



nix = tr (o-xP/s) , ruy =tT {(JyPp) , 



(2) 



where pp is the mean field canonical equilibrium state 
exp(— /3i?^]^y)/tr(exp{— /3i?^|!^y)), for rux and m^, sub- 
stitute them in H^j^j, and pp, and we are then ready to 
find the single-body physical properties of the system in 
the mean field limit. Here (3 = ^r^, with T denoting 
temperature on the absolute scale, and ks denoting the 
Boltzmann constant. 

The entanglement mean field theory begins by replac- 
ing an identity, on a site that is neighboring the inter- 
acting spins of a two-spin interaction term, by a square 
of the Pauli matrix that is involved in the interaction, 
for all the two-spin interaction terms in the Hamilto- 
nian. [An averaging needs to be done for two-spin in- 
teractions involving two Pauli matrices.] Therefore, the 
term <7^~^''cr^''' in a Hamiltonian on a two-dimensional 
square lattice can be replaced by (T^~^''CT^''cr^'+^''cr^'+^''. 
The latter can be re-written as AB, with A — cr^'''cr^'+^'', 
and B = ax~^'^<Jx'^^'^ ■ We then assume that a certain 
pair of two neighboring spins are "special" . We then re- 
place the non-special two-spin interactions (with nearby 
spins) in all the interaction terms by their mean values. 
In a c?-dimensional cubic lattice, there will be 2v terms 
(with V -\- 1 — 2d being the coordination number of the 
lattice), that will have the special pair, along with two 
other spins in two other lattice sites. The EMFT-reduced 
Hamiltonian, for the transverse quantum ferromagnetic 
XY model on a d-dimensional cubic lattice, will therefore 
be 



n 



XY 
EMFT 



= -Jv 



(l+7)C,.aX + (l-7)a 



,(3) 



where we have ignored the terms in the Hamiltonian 
which will not contribute to the EMFT equations below, 
and where we have assumed that the neighboring lattice 
sites i and j are special. Cxx and Cyy are respectively the 
XX and yy correlators, and the self-consistency equations 
(EMFT equations) are 



Cx 



tr 



913 



yy 



tr 



J 9 13 



(4) 



where gp is the canonical equilibrium state 
exp(-/?Hfl^^T)/tr(exp(-/37^f];,p^)). The EMFT 
equations are to be solved for Cxx and Cyy, and sub- 
stituted in "H^Jjprp and gp. We can then calculate the 
two-particle physical properties of the physical system 
described by H^^ in the EMFT hmit, and look for 
possible phase transitions of the system. 

A similar formalism works, with slight modifications, 
for classical spins, higher quantum spins, more complex 
lattices, etc. Also, both the mean field theory as well 
as the EMFT has been described for the ferromagnetic 



cases. The antiferromagnetic case requires some modi- 
fications in the mean field theory, and correspondingly 
some changes in the EMFT. These will not be discussed 
in this paper. 



A. The Transverse Ising model 

For simplicity, let us consider the transverse quantum 
Ising model (7 = 1). This is the simplest model which 
exhibits a quantum phase transition (at zero tempera- 
ture) 0, Q , that has been experimentally observed . 
for which the EMFT equation reads 



Cxx 



^Ising 
■'EMFT 



sinh(2/3i/C^^J) + ^^-^sinh(2/3r) ,(5) 



where the corresponding EMFT Ising partition function 
is given by Z^^^pj, = 2cos\y{2(3i'CxxJ) + 2cosh(2/3r), 

with r = ^v-^clxJ^ + K^. 

We are now ready to present the EMFT phase dia- 
gram, where we use the xx correlator, Cxx, as the order 
parameter of the many-body system. Let us first look at 
the T = picture, where transitions are driven by quan- 
tum fluctuations only. The zero temperature EMFT Cxx 
is given by 



Cxx 



1 



7^' 



(6) 



for h < he = vJ. It is vanishing otherwise. Let us now 
fix our attention on the other extreme: the behavior of 
the system with respect to temperature for zero field, in 
which case the EMFT equation reduces to 



Cxx — tanh {2l3vCxxJ) 



(7) 



whereby a temperature-driven phase transition is ob- 
tained at T = Tc = 2i'J/kB- The complete phase dia- 
gram can be seen in Fig. 1, which contains these extreme 
cases as special instances. 

It is possible to calculate the critical exponents in the 
EMFT hmit. The critical exponent for the EMFT xx 
correlator can be calculated as follows, which we find 
for both the temperature-driven phase transition on the 
h — Q axis, and for the quantum phase transition on the 
T = axis. In the zero temperature scenario, the critical 
exponent is ^, as can be found by using Eq. ([6]). In 
the zero field case, we can perform an expansion of the 
equation. 



2f3vCxxJ = tanh Cxx, 



(8) 



around T = Tc, and the critical exponent is again i. In 
a similar fashion, one can obtain the critical exponent 
for the EMFT energy gap to be i. Note that the same 
exponent is obtained for the gap in the MET limit 0]. 

The entanglement mean field theory can not only pre- 
dict the correlations of the system, but one can also study 
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FIG. 1: The xx correlator, Cxx, of the transverse Ising model 
in the EMFT limit. The temperature-driven phase transition 
at T" = Tc and the quantum phase transition at h = he are 
clearly visible on the respective axes. All the three axes rep- 
resent dimensionless quantities. 



FIG. 2: EMFT entanglement in the transverse Ising model, 
plotted against temperature and field strength. The base axes 
represent dimensionless system parameters, while the vertical 
axis represent entanglement, measured in ebits. 



the behavior of its entanglement. Apart from its fun- 
damental importance, entanglement is known to be the 
basic ingredient in quantum information tasks We 
will quantify the entanglement of a two-party quantum 
state Q, by its logarithmic negativity defined as 



E{e)^\og^ II 111, 



(9) 



where || • || denotes the trace norm of its argument, and 
Q^^ denotes the partial transpose of g with respect to one 
of the two parties forming the state g. The behavior of 
entanglement in the EMFT canonical state with respect 
to temperature and applied field, in the transverse Ising 
model, is seen in Fig. 2. 

The other members of the class of Hamiltonians given 
in Eq. ([T]) have a similar behavior with respect to their 
two-body physical properties in the EMFT limit. We 
can compare this with the fact that they fall in the same 
universality class Q ■ A different universality class is con- 
sidered in the succeeding section. 



III. THE EMFT-REDUCED XX MODEL 

The transverse field quantum XX model on a d- 
dimensional cubic lattice is described by the Hamiltonian 



H 



XX 



2 ^ 



hj2<, (10) 



where J and h are positive. The physical importance of 
the Hamiltonian includes that it can be obtained from the 
Bose Hubbard Hamiltonian for hard-core boson limit, by 
suitably associating the bosonic creation and annihilation 



operators with the Pauli matrices 5] . The corresponding 
EMFT-reduced Hamiltonian is 



n/X X 

n 



EMFT 



-JvC 



(11) 



where C is the xx (which is same as the yy) correlator of 
the system, and where we have supposed that the neigh- 
boring lattice sites i and j are special. 
The EMFT equation in this case reads 



C ^ 2sinh(2/3i/CJ) /Z^^^pT 



(12) 



where the EMFT XX partition function is given by 
Z^^pj, = 2 cosh (2/3z/C J) + 2cosh(2/3/i). Solving for C 
from the EMFT equation, we can subsequently find other 
physical quantities of the system. In particular, the en- 
tanglement of the system, as quantified by its logarithmic 
negativity, is given in Fig. 3. 

It is also possible to write an EMFT equation for 
the ground state of the EMFT-reduced XX Hamiltonian, 
solving which we find 



C = 1, for — < 1 
Jv 

C = 0, for — > 1. 
Jv 



(13) 



It is interesting to compare this result with the fact 
that the physical system represented by the Hamiltonian 
, in the one-dimensional case, undergoes a Mott in- 
sulator to superfluid transition at h/ J = 1 [sl. [Tst- 



IV. EMFT AND THE HEISENBERG MODEL: 
CORRELATION CURIE- WEISS LAW 



Investigations in magnetism in solids very often starts 
off by using the Heisenberg Hamiltonian T4] , which, for 
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FIG. 3: EMFT entanglement in the transverse XX model. 
The plot shows the behavior of entanglement of the canonical 
state of the EMFT-reduced transverse XX model with respect 
to temperature and applied field. The base axes represent 
dimensionless system parameters, while the vertical axis rep- 
resent entanglement, measured in ebits. 



a d-dimensional cubic lattice, is given by 



H 



HD 



fj* (7'^ -\- (T^ (T'^ -\- <7^ (J"' 
^ X ' y y ' ^ z 



(14) 



where J, h, and /Lt are positive. An important conclusion 
of the MFT treatment of this model is the Curie- Weiss 
Law, which predicts the behavior of magnetization of the 
physical system in its paramagnetic phase. We will see 
that it is possible to extract a similar law for the correla- 
tions in the system by solving the corresponding EMFT 
Hamiltonian. 

The EMFT-reduced Heisenberg Hamiltonian is given 

by 



The partial derivative of magnetization with respect 
to the applied field is defined as the susceptibility of the 
system, and the usual Curie- Weiss Law is given for that 
quantity, for vanishing applied field. We define the cor- 
relation susceptibility of the system as the partial deriva- 
tive of the correlation C with respect to the field strength 
h. This quantity, as solved from the EMFT equation (Eq. 
(fTB)l ). for /i = 0, is given in Fig. 4. Note that the correla- 
tion susceptibility clearly signals the onset of the param- 
agnetic phase of the system at T = Tf^ = vJ/ks- The 
data obtained for —^1, „ can be fitted to the curve 

fj. on \h—0 

(f^) a/{T - Tf^)^, for T > 1.2TfP, and it gives the 

optimal values of the curve parameters as a = 0.06895 
and A = 1.08162. The corresponding mean square error 
is 1.2 X 10-4. 
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FIG. 4: Curie- Weiss Law of EMFT-reduced correlation in the 
Heisenberg model. The partial derivative of the zz correlation 
(which is same as the xx and yy ones) with respect to the field 
h is plotted, as red triangles, against temperature, for ft = 0. 
More precisely, ^ ^ | is plotted on the vertical axis as a 

function of ^^j^ on the horizontal axis. Both axes represent 
dimensionless quantities. The blue squares represent a fit of 

the plotted data with the function (^f^^ ol/{T - Tf^)^, for 



T > 1.2T<5 



See text for further details. 



'^EMFT — 2Jl'C 



-hfi 



(15) 



where C is the zz (which is the same as the xx and yy) 
correlation of the system, and where we have supposed 
that the lattice sites i and j are special, for constructing 
the EMFT Hamiltonian. The EMFT equation in this 
case is 



C 



o-y-HD 
•^■^EMFT 



[sinhFc — exp(— Fc) -I- expFc coshF^ 



(16) 

where the EMFT Heisenberg partition function is given 



by Z_ 



HD 
EMFT 



~ 2coshrc + 2exprc coshF/i, with Fc = 



2(5vCJ and F,, = 2V3phfi. 



V. CONCLUSIONS 

The mean field theory is a useful tool to obtain im- 
portant information about single-body physical quanti- 
ties of many-body systems, especially the ones which are 
not tractable analytically. We have presented an entan- 
glement mean field theory that can be used to obtain 
information about two-body physical quantities of many- 
body systems. The theory predicts the phase diagram of 
the physical system, and the critical exponents of their 
two-body quantities, as we have shown for several im- 
portant classes of many-body systems. In particular, we 
have derived a Curie- Weiss Law for correlations in the 
Heisenberg spin model. 
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